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We analyze the effects of lateral device size and magnetic material parameters on the fer- 
romagnetic resonance (FMR) response. Results presented are directly relevant to widely used 
FMR experimental techniques for extracting magnetic parameters from thin films, the results of 
which are often assumed to carry over to corresponding nanometer-sized patterned devices. We 
show that there can be significant variation in the FMR response with device size, and that the 
extent of the variation depends on the magnetic material properties. This explains, for example, 
why different experiments along these lines have yielded different size-dependent trends from 
damping measurements. Observed trends with increasing size and different material parameters are 
explained through the evolution of three distinct eigen-modes, demonstrating the respective roles of 
demagnetization and exchange. It is also shown that there is a crossover of dominant eigen-modes 
in the response signal, accompanied by conjugating edge-type modes, leading to evident effects 
in measured linewidth and damping. Among the sizes considered, in higher saturation magneti- 
zation, we observe as much as a 40% increase in apparent damping, due solely to device size variation. 
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I. INTRODUCTION 

When a ferromagnetic material experiences a static 
magnetic field, its ground state is additionally split by 
the Zeeman energy, resulting in a condition for absorp- 
tion of electromagnetic power by ferromagnetic resonance 
(FMR) given by l^Ez = g^sBeff = hu; is the Bohr 
magneton; B^ff is the effective magnetic induction; and 
w is the frequency of the additional AC field, which is usu- 
ally in the gigahertz rangeji"— For more than 60 years, 
tapping into such a phenomenon has proven to be a pow- 
erful tool to extract information from magnetic materials, 
including both static (e.g. magnetic anisotropy, effective 
magnetization) as well as dynamic parameters (e.g. g- 
factor, damping parameter a)i^^— The damping parame- 
ter a, which quantifies the magnetization relaxation rate, 
is crucial for magnetic systems driven by finite temporal- 
length inputs, such as a current pulse or a modulated field 
on a hard-disk drive (HDD). On a HDD, a determines 
the time evolution of magnetic grains in recording media, 
and in Spin- Transfer Torque Magneto-Resistive Random 
Access Memory (STT-MRAM), it also determines the 
critical current to control magnetic states. Therefore, the 
performance of such data storage devices relies strongly 
on an accurate knowledge of damping. In FMR experi- 
ments, the concept of apparent damping is often used for 
the measured quantity to make the distinction between 
that which is measured experimentally by inferring from 
linewidth and resonance measurements and the actual 
intrinsic damping parameter a in the Landau-Lifshitz- 
Gilbert (LLG)ii equation. Apparent a is investigated 
here directly using FMR analysis. 

For some common FMR techniques, in order to obtain 



a strong enough signal in a measurement, the magnetic 
sample often needs to be orders of magnitude larger lat- 
erally than the device of interest. However, with the 
advent of STT-MRAMi^ opening the way to magnetic 
storage devices with characteristic dimensions less than 
20nm, there is a need to understand the fidelity of mea- 
surements of dynamic parameters determined by FMR 
on unpatterned thin filmsHiii. While recently developed 
techniques such as Spin- Torque FMR—"—, making use 
of microwave spin-currents or even spin-orbit coupling 
induced torquesi^, allow extracting magnetic properties 
at the device level, they are not suitable for all mag- 
netic devices (ideal if designed to flow current, and con- 
tains few magnetic layers). Other techniques, such as 
Frequency- or Time- Resolved MOKE or Magnetic Res- 
onance Force Microscopy have also been used in recent 
years to access structures smaller than typical thin films, 
which, aside from facing their own challenges'^, are also 
not as widespread. Thus, some of the FMR techniques 
in use are incapable of resolving nanometer sized device 
parameters. And, it is unclear down to what size will 
measurments remain valid. The effects of lateral size on 
the FMR response have been investigated experimentally 
in previous studies. For instance, Shaw et al. studied size 
effects in elliptical arrays and observed no change in ap- 
parent a from the patterning process^^ down to 50nm, 
while a later study by Noh et al. found that a changes 
significantly with the device size— It is worth pointing 
out that these studies use different materials. Thus, it 
is still an open question as to whether FMR measure- 
ments on larger size films yield results which represent, 
with fidelity, the device of interest, which is often much 
smaller. 
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In this work, we resolve this discrepancy through ex- 
tensive 3D finite element method (FEM) simulations, 
investigating the effects of lateral device size and mag- 
netic material parameters on the FMR response. We 
show that both results are indeed possible and demon- 
strate why this is the case. In particular, we discuss 
resonance frequency, lincwidth, and apparent a behav- 
ior with size variation. We also vary the magnetic 
material properties using the saturation magnetization, 
0.6T< HqMs < 2.4T, as it controls significant con- 
tributions to the resonance response, both through the 
demagnetization field, and equally important, through 
the exchange length l^x = \/2AJJiqM^-, A is the ex- 
change stiffness; /^o is the magnetic permeability in free 
space. The chosen range of ^^Mg covers most materials 
of interest from Permalloy {fioMg ^ IT) to CoFe alloys 
i^ioMs ^ 1.8 - 2.2T)- . As wih be shown, smaller lex 
permits higher order (i.e. more nonuniform) modal con- 
tributions into the response sooner with respect to size, 
and is important to consider. 

We first describe the FEM model used and demon- 
strate a validation of it, reproducing reported experi- 
mental and simulated data. Then, we discuss the res- 
onance behavior versus device size across different ma- 
terials. We especially emphasize the behavior of devices 
small enough for uniform precession modes to be domi- 
nant, as small as lOnm, then study the response through 
the evolution of the coexisting cigen-modes when scaling 
up size laterally and varying magnetic material proper- 
ties. In particular, our discussion of the corresponding 
resonance frequency, linewidth, and a response is given 
based on three relevant spin-wave eigen-mode evolutions. 
Finally, we summarize our results. 



II. FMR MODEL VALIDATION 

Our analysis utilizes an FEM formulation implemented 
within the software package SpinFlow 3D.— The model- 
ing procedure is as follows: first, a micromagnetic equi- 
librium state is computed with the applied DC field Hoc- 
This is followed by a computation of the spin- wave eigen- 
modes and associated eigen-frequencies in the vicinity 
of the known equilibrium state. With this information, 
the FMR response to a small AC field hrf is deter- 
mined. More details of the software formulation can be 
found clsewhcroSi. An important requirement of this ap- 
proach is to validate the model by simulating indepen- 
dently reported experimental and simulated data. For 
this purpose, we compare our results with reported data 
for Permalloy (Py) nanodisks from Shaw et al^. The 
simulated device consists of a single magnetic layer sand- 
wiched between two nonmagnetic layers, as illustrated in 
Fig. [IJa). Hue is applied along y. 

Note that the sizes indicated are nominal sizes as de- 
fined in Ref. |3 Good agreement is achieved between 
our simulation results of the resonance frequencies and 
reported experimental data for all sizes considered (Fig. 




Frequency {GHz) 



Figure 1. (Color online) (a) Device geometry used in our sim- 
ulations: a magnetic layer (dark blue) is sandwiched between 
two nonmagnetic layers (light green). For the model valida- 
tion, the geometry is elliptical, with the applied field along 
the long axis of the ellipse, as in Ref. Q. (b) A typical FMR 
spectrum showing two visible resonance peaks simulated for 
a circular disk with diameter d — 150nm, fioMs — 1.7T, and 
Hoc = 0.05T. 

[5]), confirming the validity of our model. In the following 
sections, we present our analysis based on simulation re- 
sults using the same numerical procedure to investigate 
the FMR response in more detail, varying lateral size and 
material parameters. All simulation parameters used are 
also summarized in Table I. 
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Figure 2. (Color online) Comparison of simulated resonance 
frequency response for Py devices with nominal size (a) 50nm, 
(b) lOOnm, and (c) 200nm as per Ref.fl^. Each figure includes 
experimental data (black circles), reported simulations (thin 
blue-line, NIST Simulation), and our simulation results (thick 
red- line, Our Simulation). 



III. RESULTS AND DISCUSSION 

From here out. we discuss the FMR response of cir- 
cular disk-shaped devices with varying lateral diameter 
(i, where lOnm < d < 300nm. While magnetic devices 
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Table I. Parameter values in simulations 



kNrl 



Parameter name 


Figure 2 


Other figures 


MO Ms 


0.96T 


0.6 - 2.4T 


exchange stiffness A 


llpj/m 


20pJ/m 


thickness t 


lOnm 


2nm 


mesh size 


5nm 


2nm for (f> > 40nm 






Inm for (f> < 40nm 


shape 


elliptical 


circular 


size 


a.71 X 65nm2 


10 - 300nm (diam) 




b.l20 X llOnm^ 






C.235 X 215nm2 




Gilbert damping a 


0.01 


0.01 



with features as small as lOnm are challenging to fab- 
ricate and access experimentally, numerical simulations 
avoid this difficulty and enable us to observe behavior 
and address whether thin film measurements obtain rep- 
resentative device relaxation parameters down to sufS- 
ciently small device sizes. Also note that in our FEM 
model, we do not take into consideration spin pumping 
nor inhomogeneous properties that would contribute to 
field linewidth broadening. Thus, our results reveal con- 
tributions due solely to size variation. 

A typical frequency-swept FMR spectrum xif) is 
shown in Fig. [ijb). There are, in general, two observ- 
able peaks with significant amplitude in the frequency 
range considered, indicating the presence of more than 
one eigen-modc. We consider three relevant eigen-modes 
and the resulting interplay and evolution affecting both 
observed peaks. Each resonance frequency fres is ob- 
tained by fitting each peak of the simulated FMR re- 
sponse to Lorentzians of the form 



xif) 



2P 



A./ 



7T 4(/ + /,es)2 + A/2 



(1) 



where P is the peak amplitude, / is the frequency of the 
applied AC magnetic field hrf, and A/ is the full width 
at half maximum. In Fig. [3l we plot fres versus disk 
diameter for three distinct eigen-modes which contribute 
to the response within the frequency range considered. 
Note that /ioMj = 1.7T, under an applied DC field of 
0.05T, enough to saturate the magnetization along y. In 
all illustrations of mode profiles (shown as insets in Fig. 
[3]), the color scale corresponds to the amplitude of pre- 
cession - from largest (red) to smallest (blue). From the 
spatial distribution of precession amplitude, we identify 
uniform (U), edge (E) and center (C) modes, accordingly. 
We observe that eigen-mode profiles generally evolve as 
the diameter is increased, and are labelled, for instance, 
U-E to represent evolution with size from uniform to edge 
mode. Next, we discuss the origin of the evolution of 
modes as they are intimately tied to the behavior of fres- 
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Figure 3. (Color online) Frequency response for the 
(a) Uniform-to-Edge (U-E), (b) Edge-to-Center (E-C) and 
(c) Edge-to-Edge (E-E) eigen-modes versus device size for 
^oMs = 1.7T, applied field Hoc = 0.05T. Spatial mode pro- 
files of precession amplitude are shown as insets. 



A. Evolution of Modes 

For devices with sizes on the order of lex , which we re- 
fer to as "small" devices from hereon, the lowest energy 
eigen-mode is the uniform precession mode^i. Due to 
the spatial distribution of the demagnetization field, the 
uniform-mode (U-niodc) is seen to evolve into an edge- 
mode (E-mode) that can be observed as device size is 
increased. Wc refer to this evolution as the U-E mode. 
To see the correlation more clearly, we also show the com- 
puted demagnetization field distribution in Fig.UJ^a), cor- 
responding to the magnetization equilibrium state for a 
150nm device. This distribution is typical for all sizes. 
Notice the inhomogeneity of the demagnetization field 
along the edges in the direction of Hoci which is known 
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to create a potential well for spin-waves tending to lo- 
calize them^^. The modal evolution shown in Fig. ^a.) 
is a direct result of the demagnetization energy as well 
as the distribution. To further illustrate demagnetiza- 
tion edge localization, we also plot in Fig. HJ^b) and (c) , 
respectively, computed demagnetizing field components 
parallel to the DC field (y) and Hf^ perpendicular 
to the film plane (z) for different lateral sizes. The de- 
magnetizing fields are clearly strongest at the edges, and 
comparatively weak in the interior region. Note, however, 
these observations are computed for the case of nearly 
uniform magnetization along {y). As size is increased 
sufficiently beyond lex , nonuniform affects manifest more 
strongly, leading to more complex behavior. 

A second eigen-mode evolution is shown in Fig. [3][b), 
which evolves from an edge-mode into an eigen-mode 
with dominant precession amplitude within the center 
of the device. This higher-order mode starts with edge 
precession confinement along both x (to lesser extent) 
and y ~ a four-fold symmetry. We refer to this distinct 
modal evolution as the E-C mode (edge-to-center). This 
particular mode shows the greatest nonuniformity, and 
consequently involves more exchange energy in addition 
to the demagnization field. For a magnet with nonmag- 
netic interfacial boundaries, exchange energy is known 
to be lowest at the edges, due to the free-spin boundary 
conditionf22 In this case, the system develops more inte- 
rior motion for access to more exchange energy useful to 
lower its total magnetic energy arising from demagneti- 
zation field effects. 

The third relevant eigen-mode is shown in Fig. EJc), 
in which edge-dominated precession solely persists for the 
entire frequency range considered. Therefore, it begins 
and ends as an E-mode, also localized along y edges. 
We refer to this as the E-E mode (edge- to- edge), and its 
evolutionary behavior is also explained by the demagne- 
tization field, as discussed above. Next, we analyse the 
implications of the modal evolutions on fres- 



B. Resonance Frequency Behavior 

To understand the resonance frequency response for 
non-uniform higher order eigen-modes, we consider a gen- 
eralized form of the Herring-Kittel formula^i acounting 
for quantization in all three dimensions: 



ZTT Mg Alg 

(2) 

where q is the wave vector given here by 

+(t) 

The parameters m, n, and p are quantization numbers 
arising from the finite size of the device; t is the device 
thickness; d is the lateral dimension (i.e. diameter); and 
p contains nontrivial demagnetization field informa- 
tion, which also depends on the wave vector. Note that 
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Figure 4. (Color online) (a) Demagnetization field (amplitude) 
distributions in the middle plane (along z) of the magnetic 
layer, computed from the magnetization equilibrium distribu- 
tion for fioMs = 1.7T, for a 150nm device. Color scale shows 
Il-Wijll. Computed (b) in-plane and (c) out-of-plane demag- 
netization field component along y for x=0 (center line) for 
different sizes, (d) Resulting local ferromagnetic resonance 
using Eq. (|8]). At both center and edge, the resonance fre- 
quency always increases with increasing size. 



although no simple quantization scheme can be used for 
structures confined in all three dimcnsionsii, the pre- 
sented form of the wave vector does allow us to qualita- 
tively understand the behavior of fres for higher order 
eigen-modes. 

In the size regime where the U-E mode is mostly uni- 
form, the increase in fres is understood by considering 
the decreasing (with increasing size) demagnetizing field. 
This decreasing trend can be seen in Fig. IH^b) and (c) . 
To understand why the field decreases as size increases, 
consider the exact solution for the demagnization field in 
the case of uniform magnetization, given byii 



Hd(x) 



1 

47r 



M n- 



rdS' 



(4) 



Consider a rectangular film with uniform magnetization 
M in-plane along y. In this case, only two of the six rect- 
angular boundary faces contribute to the integral in (|H), 
i.e. the two faces with normal vectors n parallel to M. 
fn the center, the value is identically due to symmetry, 
however, away from the center and only near the edges, 
nonzero values exist. Evaluating the field leads to terms 
of the following form 



Hd oc Ms^dS' 
d^ 



(5) 



By keeping the film thickness t constant, dS' scales only 
with d (instead of d x t) and the result is that the de- 
magnetization field scales as Hj) oc l/d. Thus, the de- 
magnetization field decreases with increasing lateral size. 
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The resonance frequency consequently increases because 
fres oc H^ff = — \H where both predommantly he 
along y. This accounts for the U-E mode increasing reso- 
nance frequency behavior in small device sizes. However, 
as size increases beyond ~ 50nm, the nonuniform mode 
behavior contributes more significantly. 

In this case, equation ^ suggests exchange energy 
also contributes to fres in higher-order modes, through 
the wave- vector q. In particular, it follows from ([2]) that 
resonance frequencies generally decrease with increasing 
device size, as \/d?. This is in qualitative agreement 
with trends shown in Figures [31[b) and (c) corresponding 
to higher order eigen-modes. This also applies to Fig- 
ures [3{a), after the transition to an edge mode has taken 
place. In the nonuniform modes, the gradient of the 
eigen-modes also contains exchange information as they 
evolve with size. The modes indicate precession ampli- 
tude, propertional to normal components of the magne- 
tization (normal to y). The total exchange energy given 
by 

^A = W^j i^^^fdV (6) 

By observing the trend in the gradient of the mode, the 
trend of the exchange energy can also be seen. Fig- 
ure [5]^a) illustrates the qualitative trend of the exchange 
associated with the E-C mode (the most nonuniform 
mode), indicating on the eigen-mode vector directions 
(with black arrows). As size increases, the complexity of 
the mode reduces, and this generally reduces exchange 
energy density. This behavior is equivalent to reducing 
the quantization numbers m and n. For two of the device 
sizes in the region of sharp decrease in f^es , the frequen- 
cies are plotted in Figure [5Kb), along with the gradient 
distributions of the in-plane normal mode component (x) , 
and its maximum computed gradient indicated (in arbi- 
tray units). This demonstrates an overall decrease in the 
exchange energy through the E-C mode evolution, and 
this contributes to a reduction in fres for higher order 
modes. 
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Figure 5. (Color online) (a) Illustration of reduction in ex- 
change energy in high frequency E-C mode, with increasing 
size (50, 150, and 300nm shown) and (b) Computed gradients 
of in-plane normal modal component for 50nm and lOOnm. 
Maximum gradient (in arbitrary units) shown beside the re- 
spective gradient distribution. 



for the dominant mode of the FMR response signal as 
size increases. The cross-over is illustrated in Fig. |6j 
seen clearly for fXoMs = 1.7T and 2.4T. 

Near the cross-over size, there also occurs a conjuga- 
tion or joining of the U-E and E-E eigen-modes, as shown 
in Fig. [71 where fi^Mg = 1.7T. At the point of conjuga- 
tion, the mode profiles of the U-E and E-E modes be- 
come similar or quasi-degenerate and resonate at similar 
frequencies. Both the cross-over and the conjugation of 
modes can have effects on the measured linewidth as will 
be discussed in a later section. 



C. Cross-over and Conjugation of Modes 

Another important characteristic of the FMR spec- 
trum, aside from the resonance frequency trend with size, 
is the relative amplitudes of the observed modes. The 
amplitude of an FMR peak is proportional to the volume 
of the device participating in the precession. For exam- 
ple, in the U-E mode, the edge precession is increasingly 
confined to the edges with increasing size (as can be seen 
in the insets of Fig. [3]) and the relative amplitude of the 
peak decreases, with additional contributions from a re- 
ducing demagnetization field. In contrast, the amplitude 
of the E-C mode increases with size relative to the other 
modes. Therefore, in smaller sizes, the U-E mode tends 
to dominate the signal, however, as size increases, the 
U-C mode begins to dominate. This leads to a cross-over 



D. Effect of Different 

One key observation in this study is that the FMR 
response can be sensitive to magnetic material parame- 
ters. For example, the effect of different Mg on fres can 
be seen in Fig. [8] For the U-E mode, the device size at 
which an initially increasing fres transitions to decreas- 
ing fres (due to the onset of the edge mode) is smaller in 
materials with larger Ms- For the smallest ^qMs (0.6T), 
a decrease in fres is not observed at all in our range of 
sizes. However, devices with higher have shorter lex, 
and it follows that non-uniform modes may contribute 
more significantly in smaller lateral sizes. Hence, the ef- 
fect of exchange, which tends to lower fres as discussed 
above, appears earlier. Comcomitent with this, we also 
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Figure 8. (Color online) (top) Frequency response of three 
observed eigen-modes for (a) fioM^ = 0.6T, (b) LOT, (c) 1.7T 
and (d) 2.4T. 



Figure 6. (Color online) (a) Simulated FMR response versus 
device size for ^oMs = 0.6T, (b) LOT, (c) 1.7T and (d) 2.4T. 
White arrow indicates the location of the change in domi- 
nant modes, crossing over from U-E mode to E-C mode with 
increasing size. 
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Figure 7. (Color online) Resonance frequency versus size for 
the U-E and E-E modes for noMs = 1.7T, zooming in to 
the point of conjugation. Note the increasingly similar mode 
profiles of the U-E and E-E modes. 



observe the occurence of both the modal conjugation and 
dominant mode cross-over at smaller sizes as Ms is in- 
creased. The effects on the cross-over is evident in Fig. 
M.c) and (d). 



E. Linewidth and Damping Behavior 

The linewidth A/ for each size is obtained by tak- 
ing the full width at half maximum of each Lorentzian 
peak. For a multi-peak spectrum, we can identify sepa- 



rate linewidths for each peak. We typically see two sig- 
nificant peaks which we identify as "low" and "high" fre- 
quency peaks. The linewidth versus size behavior for the 
high frequency peak is mainly determined by the E-C 
mode. The linewidth for the low frequency peak can be 
more complex because both U-E and E-E modes can con- 
tribute to the signal due to their conjugation. However, 
below the conjugation point, both the resonance frequen- 
cies and relative amplitudes of the U-E and E-E modes 
are very dissimilar, and the low frequency peak linewidth 
is controlled mainly by U-E. As size increases, the rela- 
tive amplitude of the U-E mode decreases and becomes 
comparable to the amplitude of the E-E mode. In some 
parts of this regime, their resonance frequencies do not 
coincide exactly. When the U-E and E-E fres modes do 
not exactly match, but are close enough, there is some ar- 
tifical broadening of the measured linewidth. Once their 
resonance frequencies coincide, this artificial broadening 
is lost. To illustrate this, we plot the linewidth of the 
low frequency peak (which includes contributions from 
both U-E and E-E modes, corresponding more closely to 
a measured signal) and the linewidth of only the U-E 
mode in Fig. [HI The insets show spectra of the U-E and 
E-E modes in the size range where the transient broad- 
ening may occur. Note the increased linewidth of the 
low frequency peak at 200nm which can be attributed to 
modal conjugation. 

For experiments, however, looking at the dominant 
peak (i.e. with the highest intensity) may be more rele- 
vant because limitations in detection sensitivity may not 
be enough to resolve peaks other than the dominant one. 
Without the aid of rigorous calculation or simulation, it 
could be difficult to trace the origin of a particular peak 
in the signal; our results show that the dominant peak 
may be a different mode for different sizes. In the fol- 
lowing analysis of linewidth behavior versus size for the 
dominant peak, we use the field linewidth A.H, which is 
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Figure 9. (Color online) Linewidth versus size for the low 
frequency peak (blue squares: with contributions from both 
U-E and E-E modes) compared to the linewidth of only the 
U-E mode (black circles). Broadening from conjugation is 
seen at 200nm. Inset shows the FMR spectrum for U-E (thin 
black) and E-E (thick red) modes in the size range where 
broadening may occur. 



obtained from A/ via 



A/ 



(Sf/SH) 



(7) 



The denominator is calculated from the general Kittel 
equation^ 



f ^ {H + Hi){H + H2 



(8) 



where H is the applied DC field and Hi and H2 are 
stiffness fields. We plot the linewidth of the dominant 
peak for different Ms values in Fig.[TUl We observe a clear 
increasing trend, with slightly more complex behavior for 
HqMs = 2.4T. Because the linewidth is given by 



AH = Ai7n 



7 



(9) 



and is therefore proportional to /res, we find that the 
change in dominant mode, which is accompanied by a dis- 
tinct jump to higher fres, results in an abrupt increase in 
AH at around 200nm. Moreover, after the cross-over, if 
the linewidth follows /res, the now-dominant E-C mode 
resonance frequency decreases, so the linewidth should 
decrease according to Eq. |9l where AHq is the inhomo- 
geneous linewidth broadening due to inhomogeneities in 
the device parameters, e.g. dead layers and edge dam- 
age, which are not considered here. However, linewidth 
does not decrease with device size, as does fres- This 
indicates that a is non-constant and is contributing to 
the observed linewidth behavior. 

We therefore extract the apparent damping a for each 
peak - by fitting applied field-dependent linewidth data 
to Eq. ini - to eliminate contributions from fres ■ The 
results are shown in Fig. [TT] for different Ms- For the 
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Figure 10. (Color online) Linewidths from FMR response 
versus device size for different Ms- 



low frequency peak, the increase in a is more significant 
in larger Mg. For the high frequency peak, a decreases 
monotonically for all Afg, approaching a = 0.01. 

To further explain , we consider that each differen- 
tial volume in the device experiences a slightly differ- 
ent Hef f and therefore has a slightly different fres com- 
pared to its neighbor. This difference is more pronounced 
as Ms increases because shorter lex allows more non- 
uniformity and the effect of the demagnetization field 
is also stronger. The collective response - the sum of 
these individual responses which do not necessarily over- 
lap - leads to broading. The low frequecy peak (mainly 
from the U-E mode) goes from uniform to non-uniform, 
leading to more broadening with size. However, for the 
the higher frequency peak (E-C mode), which becomes 
more uniform with increasing size, broadening generally 
decreases. 

The behavior of the linewidth for /ioA/s = 2.4T in 
larger sizes is actually explained by the corresponding 
decrease in a at 200nm as dominance shifts from the low 
frequency peak to the high frequency peak. This also 
contributes to the subsequent increase at 300nm, as can 
be seen in Fig. [10] 

The effects of size and material properties can also be 
seen clearly in Fig. [TTJ In the range considered, mate- 
rials with low Ms generally show neglible variation in 
damping, while for higher values of Ms, there is a clear 
variation with size, indicating that the extent of the vari- 
ation is sensitive to material properties, which in our 
case is Ms- In the present case, the characteristic length 
that controls the FMR behavior with size is lex- Thus, 
our results suggest that the different findings on a ver- 
sus size behavior-i^ii^ are not inconsistent, but stem pri- 
marily from the use of different materials. Our results 
arc consistent, quanlitatively, with both reports, where 
Py {fj.oMs = LOT) shows little variation, however CoFe 
{fioMs = 1.7T) does shows a visible increase with size. 

Therefore, since apparent a may be non-constant with 
size, care must be taken in using film measurements for 
device-related parameters. For a multi-peak FMR re- 



8 



0.0150 



0.0125 - 



0.0100 



0.0150 



0.0125 - 



0.0100 



n M (T) Low Frequency Peak 



^0 s 

0.6 
1.0 
-A- 1.7 
■▼- 2.4 



■ = : = = :§ 





(T) 




1.0 




1.7 




2.4 



High Frequency Peak 



50 100 150 200 250 300 
Device Size (nm) 

Figure 11. (Color online) Apparent a from FMR response for 
different Ms- Dashed line in (a) for fioAIs — 2 AT is a guide 
to the eye, highlighting the increased non-linearity around 
200nm. 



sponse, a of the dominant peaks in the extremes of the 
size spectrum considered here seem to converge to the 



same value, however, they do so on opposite ends of the 
size range. Take for example, 40nm where the low fre- 
quency peak is dominant, and 300nm where high fre- 
quency peak is dominant. Both values approach a = 0.01 
for all Ms values considered. Although the applicability 
of film measurements to the device level seems possible, 
this, of course, assumes proper accounting of all other 
contributions to the linewidth. So, there can also be close 
agreement between damping measurements from film and 
from devices below a certain size in some materialsi^ 



IV. CONCLUSION 

We have investigated the effects of device size and 
magnetic material parameters on the FMR response us- 
ing a powerful 3D FEM tool. We identified three dis- 
tinct eigen-modes, the profiles of which evolve versus size 
and consequently lead to evident effects on the FMR re- 
sponse. Two of the three modes (U-E and E-E) were 
shown to conjoin at a given size and concomitant with 
this conjugation, a cross-over in terms of mode domi- 
nance is also seen. These phenomena may affect mea- 
sured linewidth and damping. The amount of sensitivity 
to size was shown to increase as Ms is increased and may 
be the reason behind the discrepancy from different in- 
vestigations, which used different materials. Within the 
size range of lOnm to SOOnm, we have observed as much 
as a 40% increase in apparent a, due purely to size ef- 
fects. Our results suggests that an FMR measurement 
can indeed lead to infiation of the intrinsic damping pa- 
rameter, depending on the material used. 
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